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A method of processing geophysical data 

The present invention relates to a method of processing geophysical data, in particular 
to a method of processing geophysical data to determine information relating to 
azimuthal anisotropy in the geophysical data. Information about the earth's interior, for 
example relating to fractures within the earth or reflector dip, may then be obtained 

Figures 1 and 2 illustrate the occurrence of azimuthal anisotropy in seismic data, one 
type of geophysical data. Figure 1 illustrates a seismic surveying arrangement in which 
seismic energy emitted by a seismic source S is received at a seismic receiver R, for 
example a geophone. The seismic source S and the seismic receiver R are constrained 
to lie diametrically opposite one another on the circumference of a circle C. The angle 
^is the angle between the source-receiver direction (on the earth's surface) and a 
predetermined direction (north in Figure 1), and is lcnown as the "azimuthal angle" or 
just "azimuth". 

In this survey, seismic data are obtained for an initial value of the azimuth. The 
source S and the receiver R are then moved around the circle C, while still be being kept 
diametrically opposite, so that the azimuth is changed to anew value, Further 
seismic data are then acquired at the new azimuth jft. This process may be repeated, 
and seismic data acquired at further azimuths fa (fc etc. 

If the elastic properties of the earth's interior were isotropic and laterally homogenous, 
the seismic data acquired by the seismic surveying arrangement of Figure 1 should be 
independent of the azimuth ^. It is commonly found however, that the seismic data 
depend on the azimuth - that is, the seismic data are azimuthally anisotropic. This 
indicates that the elastic properties of the earth's interior are not isotropic and/or are not 
laterally homogenous. Azimuthal anisotropy of the earth's interior may arise, for 
example, as a result of vertically aligned fractures within the earth. 
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In Figure 1 the earth's interior is schematically represented as an upper, isotropic layer 1 
and a lower, azimuthally anisotropic layer 3, separated by a horizontal interface 2. The 
lower layer 3 is azimuthally anisotropic owing to the presence of vertically-aligned 
fractures 4. 

Figures 2(a) to 2(c) show a schematic representation of seismic data acquired by the 
seismic surveying arrangement of Figure 1. Figure 2(a) shows the vertical component 
of the seismic energy acquired at the receiver R, Figure 2(b) shows the transverse 
component of the seismic energy acquired at the receiver R, and Figure 2(c) shows the 
radial component of the seismic energy acquired at the receiver R. (The radial 
component of seismic energy acquired at the receiver R in Figure 1 is the horizontal 
component in the source-receiver direction. The transverse component of seismic 
energy is perpendicular to the radial component and lies in the horizontal plane.) 

Each of Figures 2(a) to 2(c) shows a set of seismic energy traces acquired at different 
source-receiver azimuths. Each set contains 36 traces, and the azimuth at which one 
trace in a set is acquired is 10° greater than the azimuth at which the previous trace in 
that set was acquired. Each trace was acquired by actuating the source S to emit a 
single pulse of seismic energy, and each trace shows the amplitude of the respective 
component of the seismic energy received at the receiver Rasa function of time after 
the source S was actuated. The first event in the traces of Figures 2(a) and 2(c) 
corresponds to the Pp-wave reflection and the second event corresponds to interfering 
arrivals of the Psl- and Ps2-wave reflections. The traces in Figure 2(b) show only an 
event corresponding to interfering arrivals of the Psl- and Ps2-wave reflections - the Pp 
reflection does not appear in the transverse component. The seismic data acquired by 
the receiver R would in practice contain a "direct event" corresponding to seismic 
energy that was transmitted from the source S direct to the receiver R through the upper 
layer 1, but this has been omitted from Figures 2(a) to 2(c) for clarity. 

It will be noted from Figure 2(b) that the transverse component of the seismic data 
exhibits azimuthal anisotropy. Azimuthal variations are clearly evident in both the 
magnitude and the polarity of the event in the transverse component. The vertical and 
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radial components of the seismic data are, however, substantially independent of 
azimuth. 

In general, the following attributes of acquired seismic data can show dependence on 
azimuth. 

1) Arrival time. If the path of seismic energy from the source S to the receiver R 
passes through an azimuthally anisotropic medium then the propagation velocity 
within that medium will be azimuth-dependent. This leads to azimuth-dependent 
arrival times. 

2) Polarisation angle. In the case of the anisotropic layer 3 of Figure 1 having 
vertically aligned fractures, symmetry planes will exist in the directions normal to 
the fracture and parallel to the fracture (the "fracture strike" direction). For events 
derived from seismic energy propagating parallel to one of these symmetry planes, 
the polarisation of the events must be aligned either parallel or perpendicular to this 
plane. For seismic energy propagating in other directions, however, the polarisation 
is not similarly constrained. Thus, the deviation of the Ps-wave polarisation from 
the source-receiver azimuth would be zero in the fracture normal and fracture strike 
directions but would generally be non-zero for other azimuths. 

3) Transverse energy. The magnitude of the transverse component of seismic energy 
recorded at the receiver R in the seismic surveying arrangement of Figure 1 is 
related to the symmetry planes. For events derived from seismic energy propagating 
parallel to one of the symmetry planes no energy will be observed in the transverse 
energy component, whereas events derived from seismic energy propagating in 
other directions will have a non-zero transverse component. Moreover, in some 
survey geometries transmitted P-wave arrivals can be examined in a similar way - 
for example if direct transmitted arrivals are being recorded in a Vertical Seismic 
Profile (VSP) survey, then the P-wave arrivals can also have a transverse 
component. 
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4) Polarity - the polarity on the transverse component of the Ps events will reverse as 
the source-receiver azimuth passes through the fracture normal and fracture strike 
directions. If transmitted arrivals are being recorded, for example in a VSP survey, 
then polarity changes may also be observed in the P-wave arrivals, as reported by 
GJRossi, A.Vesnaver and S.A.Petersen in "Anisotropy Detection in a 3D three- 
component VSP", First Break Vol.19, No. 4, ppl91-200 (2001). 

5) Temporal frequency. Azimuthal anisotropy can lead to two S-events being excited, 
and these are known as the Psl and Ps2 events. Where the time delay between the 
arrivals at the receiver of these two s-events is insufficient to separate the S-waves 
into two distinct arrivals, interference between the two shear wave arrivals may 
occur. Depending on whether it is constructive or destructive, this interference will 
lead to either an increase or a decrease in the temporal frequency. Since the shear- 
wave splitting is azimuthally dependent then it follows that the temporal frequency 
is also azimuthally dependent. Azimuthal variations will also be apparent in 
temporal frequency because of azimuth-dependent attenuation. 

5) Amplitude. The amplitude of seismic energy recorded at the receiver will vary as a 
function of azimuth because of azimuthally anisotropic attenuation and wavefront 
divergence. 

6) Linearity. For the Ps events, maximum linearity is observed for events derived from 
seismic energy propagating parallel to one of the symmetry planes because only one 
s-wave will be excited. The linearity will be lower for events derived from seismic 
energy propagating in directions other than the symmetry planes. 

As noted above, azimuthal anisotropy is often attributed to the presence of vertically 
aligned fractures in the earth. Thus, characterising azimuthal anisotropy in seismic data 
or other geophysical data can provide a better understanding of the presence of fractures 
within the earth, and this can lead to improved productivity from hydrocarbon and 
geothermal reservoirs. 
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A further cause of azimuthal anisotropy in geophysical data is that many geological 
reflectors within the earth are not horizontal but are dipped. Where geophysical data is 
acquired from a dipped reflector, some attributes of the data will show azimuth- 
dependent variations. These azimuth-dependent variations will possess mirror 
symmetry planes aligned with the reflector dip and the reflector strike. Thus, 
determination of the mirror symmetry planes in the geophysical data identifies the 
direction in which the dipped reflector is oriented. 

Previously, azimuthal anisotropy in geophysical data has been investigated by analysing 
shear-wave splitting or fitting elliptical variations to normal moveout (NMO) velocities 
or AVO as a function of azimuth. 

Ellipse fitting has previously been applied to the analysis of azimuthal variations in 
seismic data attributes, as reported in US patent No 5 508 973 and by V.Grechka and 
LTsvankin in "3D description of normal moveout in anisotropic media", SEG Expanded 
abstracts ppl487-1490 (1996). R.Bale et al. have proposed, in "Analysis and stacking 
of 3D converted wave data in the presence of azimuthal anisotropy**, SEG Expanded 
Abstracts pp 1189-1192 (2000), a method of identifying polarity changes which are 
expected to occur every 90 degrees of azimuth on the transverse component for Ps 
waves. This technique applies a numerical search using a polarity filter to find an 
azimuthal angle where the polarity reversals occur. Garotta and Granger have proposed, 
in "Acquisition and processing of 3Cx3D data using converted 3Cx3D data using 
converted waves", SEG, pp995-997 (1988), a similar method that uses measurements of 
the ratio of transverse to radial components to identify 'dimming' directions. These 
dimming directions are associated with the fracture normal and fracture strike azimuths 
that occur every 90 degrees. The theoretical basis for these prior art ellipse fitting 
methods assumes that higher-order terms can be neglected. In real geophysical data, 
however, higher-order terms may be significant. 

A first aspect of the present invention provides a method of processing geophysical 
data, the data comprising at least a first set of geophysical data acquired with a first 
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source-receiver azimuth and a second set of geophysical data acquired with a second 
source-receiver azimuth different from the first source-receiver azimuth, the method 
comprising; (a) selecting an attribute of the geophysical data that is dependent on the 
source-receiver azimuth; and (b) determining a value of the source-receiver azimuth 
about which the selected attribute exhibits mirror symmetry. This provides an improved 
method of determining the azimuths of the mirror symmetry planes within the earth. 

The method may be carried out on pre-existing geophysical data. Alternatively, the 
method may comprise an initial step of acquiring- at least a first set of geophysical data 
with a first source-receiver azimuth and a geophysical set of seismic data with a second 
source-receiver azimuth different from the first source-receiver azimuth. 

A second aspect of the invention provides a method of determining an azimuthal 
symmetry plane within the earth from geophysical data comprising at least a first set of 
geophysical data acquired with a first source-receiver azimuth and a second set of 
geophysical data acquired with a second source-receiver azimuth different from the first 
source-receiver azimuth, the method comprising the steps of: a) selecting an attribute of 
the geophysical data that is dependent on the source-receiver azimuth; and determining 
a value of the source-receiver azimuth about which the selected attribute exhibits mirror 
symmetry thereby to locate a mirror symmetry plane within the earth. 

The geophysical data may be seismic data. 

A third aspect of the invention provides an apparatus for processing geophysical data 
comprising at least a first set of geophysical data acquired with a first source-receiver 
azimuth and a second set of geophysical data acquired with a second source-receiver 
azimuth different from the first source-receiver azimuth, the apparatus comprising; 
means for selecting an attribute of the geophysical data that is dependent on the source- 
receiver azimuth; and means for determining a value of the source-receiver azimuth 
about which the selected attribute exhibits mirror symmetry. 
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A fourth aspect of the invention provides an apparatus for determining an azimuthal 
symmetry plane within the earth from geophysical data comprising at least a first set of 
geophysical data acquired with a first source-receiver azimuth and a second set of 
geophysical data acquired with a second source-receiver azimuth different from the first 
source-receiver azimuth, the apparatus comprising: a) means for selecting an attribute of 
the geophysical data that is dependent on the source-receiver azimuth; and b) means for 
determining a value of the source-receiver azimuth about which the selected attribute 
exhibits mirror symmetry thereby to locate a mirror symmetry plane within the earth. 

In a preferred embodiment the apparatus comprises a programmable data processor. 

The present invention also provides a storage medium comprising a program for a data 
processor of an apparatus as defined above. 

Preferred embodiments of the present invention will now be described by way of 
illustrative example, with reference to the accompanying Figures in which: 

Figure 1 is a schematic view of a seismic surveying arrangement in which the acquired 
data may exhibit azimuthal anisotropy; 

Figure 2(a) shows the vertical component of the seismic energy acquired at the receiver 
in the seismic surveying arrangement of Figure I; 

Figure 2(b) shows the transverse component of the seismic energy acquired at the 
receiver in the seismic surveying arrangement of Figure 1; 

Figure 2(c) shows the radial component of the seismic energy acquired at the receiver in 
the seismic surveying arrangement of Figure 1 ; 

Figure 3(a) illustrates the azimuthal dependence of a seismic attribute acquired at the 
receiver in the seismic surveying arrangement of Figure 1 ; 
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Figure 3(b) illustrates the effect of a rotation of 90° on the seismic energy pattern of 
Figure 3(a); 

Figures 4(a) and 4(b) are schematic illustrations of the principle of a method according 
to the present invention; 

Figure 5 is a schematic illustration of the principal steps of a method according to the 
present invention; 

Figure 6(a) is a schematic plan view of a seismic surveying arrangement suitable for use 
with the present invention; 

Figure 6(b) is a schematic plan view of another seismic surveying arrangement suitable 
for use with the present invention; and 

Figure 7 is a block schematic diagram of an apparatus according to the present 
invention. 

The principal steps of a method according to an embodiment of the present invention 
are schematically illustrated in Figure 5. In this embodiment the geophysical data are 
seismic data. 

Initially, at step 1, seismic data are acquired. The seismic data contains sets of data 
acquired at different source-receiver azimuths. The seismic data acquired in step 1 is 
shown in Figure 5 as comprising a set of seismic data for every 10° increment in the 
source-receiver azimuth, but the sets of seismic data could be acquired using smaller or 
greater increments in the source-receiver azimuths. In the absence of any a prior 
knowledge about the azimuthal orientation of symmetry planes within the earth, the 
geophysical data preferably contains dense azimuthal sampling over the full 360° 
azimuth range. However, if there is pre-existing information about the orientation of 
symmetry planes it may be sufficient to acquire data having a smaller azimuth range 
that includes the expected azimuth of a symmetry plane. 



WO 03/009005 



PCT/GB02/03354 



9 

It should be noted that the invention can alternatively be applied to pre-existing seismic 
data having a suitable azimuth coverage. In this case, the step of acquiring seismic data 
is not required. 

Next, at step 2, the seismic data is processed to determine one or more attributes of the 
data that are dependent on the source-receiver azimuth. One attribute that is found to be 
azimuth-dependent is selected, and the value of this attribute determined for all 
available azimuthal angles - i.e., the value of this attribute in each seismic data set is 
determined. The attribute can be any attribute of the seismic datalhat is dependent on 
the source-receiver azimuth, and could be, for example, one of the attributes (1) to (7) 
listed above. 

In general the acquired seismic data will contain data only for certain azimuths, for 
example for 10° increments of azimuth in the case of the data of Figures 2(a) to 2(c). 
Thus, analysis of the seismic data will produce a set of discrete values of the selected 
attribute. 

The results of processing the seismic data at step 2 are stored, and are thai output as 
step 3. The results may be, for example, output on a suitable display for inspection by 
an operator. It will be noted in Figure 5 that the value of the selected attribute has been 
interpolated between the discrete values derived from the seismic data to generate a 
continuous trace for the attribute over the 0° to 360° azimuth range, but this does not 
need to be done. 

In principle, the processing step described below could be carried out on the stored data 
without displaying the stored data Step 3 may therefore be omitted if desired. 

At step 4, the azimuth-dependence of the selected attribute is processed to determine 
one or more azimuthal angles about which the attribute is mirror-symmetric. The or 
each such azimuthal angle corresponds to a plane of symmetry within the earths 
interior. 
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As noted above, aligned vertical fractures within the earth's interior are one cause of 
symmetry planes within the earth. A vertical fracture causes two symmetry planes. 
One symmetry plane is in the fracture strike direction, and the other is perpendicular to 
the plane of the fracture - so that the two symmetry planes are perpendicular to one 
another. Vertical fractures within the earth are typically caused by stresses which are 
approximately constant over a large area. As a result, in a set of fractures within the 
earth the fractures will generally have similar orientations to one another, so that the 
symmetry planes of different fractures will also be similar to one another. 

The processing of the azimuth-dependence of the attribute at step 4 can be performed 
using any suitable mirror symmetry-detection technique. There are many 
known methods for mirror symmetry detection, and mirror symmetry detection is a 
well-established technique in other disciplines such as, for example, structural 
chemistry or image processing. 

Finally, at step 5, the results of step 4 are output. In Figure 5 the results are shown as 
being output as a visual trace, but the results can be output in any convenient form. 
Alternatively, the results may be stored for use in further processing of the seismic data, 
in which case step 5 may be omitted. 

Figure 5 illustrates the results of a preferred method of determining azimuthal angles at 
which the azimuth-dependence of the attribute exhibits mirror symmetry. This method 
is to generate a function indicative of the variation between the value of the selected 
attribute at a first source-receiver azimuth and the value of the selected attribute at a 
second source-receiver azimuth obtained by reflecting the first source-receiver azimuth 
about a trial azimuth of a mirror symmetry plane. The values of the attribute at these 
two azimuths are combined into a function known as an "objective function". The 
absolute value of the objective function will be at a minimum, possibly zero, if the trial 
azimuth of the mirror symmetry plane is correct. Thus, the true azimuth of the mirror 
symmetry plane can be found by finding a value of the trial azimuth of the mirror 
symmetry plane that minimises the absolute value of the objective function. This can be 
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done, for example, by finding a value of the trial azimuth of the mirror symmetry plane 
that makes the absolute value of the objective function zero. Alternatively, it may be 
done by finding a value of the trial azimuth of the mirror symmetry plane that makes the 
absolute value of the objective function less than a pre-set threshold. In this 
embodiment, the data output at step 5 of Figure 5 could be the value of the objective 
function over a pre-selected azimuth range (as shown in Figure 5), or alternatively or 
additionally they could be an azimuth of the mirror symmetry plane that minimises the 
absolute value of the objective function. 

The objective function output at step 5 of Figure 5 has a value of zero at (j> = 60°. As 
will be explained below, this indicates that the azimuth of one symmetry plane is 60°. 
Moreover, since the two symmetry planes are orthogonal to one another, it follows that 
the azimuth of the second symmetry plane is 150°. (The minimum at 0= 15° in the 
objective function arises owing to noise in the seismic data and does not represent a 
symmetry plane.) 

Determination of the symmetry planes provides information about the orientation of 
fractures within the earth, if these are responsible for the azimuthal anisotropy. 
Alternatively if the azimuthal anisotropy is a consequence of reflector dip, then the 
determination of the symmetry planes provides information about the reflector dip and 
reflector strike, and this information can be used to correct for the reflector dip during 
further processing of the seismic data. 

A preferred method of determining the azimuthal angles at which the azimuth- 
dependence of the selected attribute of the seismic data exhibits mirror symmetry will 
now be described. This method will be described using the transverse component of the 
energy acquired at the receiver as the azimuth-dependent attribute, but the method is not 
limited to this attribute. 

Figure 3(a) shows a typical result for the magnitude of the transverse component of the 
received energy, plotted against the source-receiver azimuth <f>. The magnitude is 
plotted as a polar plot, and the radius at a given azimuth denotes the magnitude of the 
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transverse component of the energy at that azimuth. The transverse component of the 
received energy displays minima at 0= 0°, 90°, 180° and 270°, and displays a 
maximum between each two adjacent minima. The magnitude of the transverse 
component of the energy plotted in Figure 3 has the same value at the points 1, 2, 3 and 
4, so that the radius of the attribute variation curve is equal at the azimuths 
corresponding to these points. The transverse energy attribute in azimuthally 
anisotropic seismic data will generally demonstrate behaviour similar to that shown in 
Figure 3(a). 

Figure 3(b) shows, in dotted lines, the result of rotating the plot of Figure 3(a) by 90°. 
This is overlaid on the original plot (shown in full lines) for ease of comparison. It can 
clearly be seen that the transverse component of the received energy shown in Figure 
3(a) does not have 4-fold rotational symmetry, since rotating the plot by 90° does not 
transform it back onto itself. Rather, the transverse component of the received energy 
shown in Figure 3(a) has two mirror planes of symmetry, and these are aligned along A- 
A' and B-B\ The mirror planes of symmetry are orthogonal to one another. 

In the case of an attribute/that exhibits azimuthal anisotropy, the existence of two 
orthogonal mirror symmetry planes requires that the following equation holds: 

/^-^o) = /(^ + <*o) = /a80-^ + ^o) = /(180 + ^-^ 0 ) (1) 

where /($is the value of the attribute at an azimuth fa and one symmetry plane lies at 
an azimuth It can be seen that this equation is satisfied for the results of Figure 3(a). 
For example, the mirror symmetry plane B-B' has an azimuth = 0°, in which case 
equation (1) requires lhat/tyf) =/0 80 + <0) = /{l 80 - f$) - and this is true for the 
results of Figure 3a. 

Where the value of is unknown, its value may be found by determining the value of 
that minimises an objective function such as the following: 
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£=Zcr i -/[-^+^]) 2 +(/;-/ii8o- ( *+^]) 2 +c/;-/[i8o+^-^]) 1 ( 2 ) 

w 

In equation (2),fi is the i * value of the selected attribute, which lies at an azimuth angle 
fa and/[a] is the value of the attribute at the azimuth a. The objective function E will 
have minimum at values of fa that are the azimuth of a mirror symmetry plane. The 
objective function will have a minimum at the azimuth of each symmetry plane, and so 
will have minima at four azimuths in the 0 to 360° azimuth range in the case of two 
orthogonal mirror symmetry planes. 

It will be noted that the objective function of equation (2) cannot be negative. This 
feature is preferred, since it simplifies use of an iterative technique to determine the 
value of fa that minimises the objective function. 

Figures 4(a) and 4(b) are a schematic illustration of the operation of this method. In 
Figure 4(a), the azimuth fa has been chosen as the first estimate of the azimuth of a 
mirror symmetry plane. This corresponds to the orientation of the mirror symmetry 
planes shown by broken lines. The open circle denotes the pointy , which is the value 
of the attribute at an azimuth angle fa. 

The value of the objective function E for the trial estimate fa of the azimuth of a mirror 
symmetry plane is found by computing the value of the attribute / at the three azimuths 
derived from the trial estimate fa of the azimuth of the mirror symmetry plane and the 
azimuth fa of the pointy. The three azimuths give the three points shown as crossed 
circles in Figure 4(a). The objective function E is then calculated from f t and from the 
value of the attribute at the three points shown as crossed circles. 

It should be noted that in a practical case in which seismic data was acquired for 
increments o£ say, 10° or 15°, then for an arbitrary initial values of fa and fa, it may 
well be found that the acquired seismic data does not contain a data set for one or more 
of the azimuths (-fa + fa), (180 + fa) and (180 + fa - fa). In this case it is necessary 
to obtain the values of the attribute at one or more of these azimuths by interpolation. 



WO 03/009005 



PCT/GB02/03354 



14 

In principle it is possible to interpolate the values of the attribute derived from the 
seismic data at step 2 of Figure 5 to generate a continuous or quasi-continuous set of 
values for the attribute, so that the output at step 3 is a continuous or quasi-continuous 
set of values for the attribute. In practice, however, it may be more efficient to perform 
the interpolation only as and when required during the calculation of the objective 
function so as to reduce the computing power required 

As shown by equation (1), the value of the attribute at the azimuth denoted by the open 
circle (# and the value of the attribute at each of the three azimuths denoted by the 
crossed circles should be equal, if the trial estimate, fa, of the azimuth of a mirror 
symmetry plane is correct This would lead to a value of zero for the objective function 
E. This is clearly not the case in Figure 4(a), however, and this indicates that the trial 
estimate fa of the azimuth of a mirror symmetry plane is not correct. The value of the 
objective function E for the case of Figure 4(a) will not be a minimum. 

Figure 4(b) corresponds to Figure 4(a), but is for a different trial estimate, fa, of the 
azimuth of a mirror symmetry plane. The estimate fa is much closer to the correct 
azimuth, and it is clear that the value of the attribute at the azimuth denoted by the open 
circles and the value of the attribute at each of the azimuths denoted by the crossed 
circles will be approximately equal to one another, so that the value of the objective 
function E will be close to a minimum. 

One way of determining the azimuth of the mirror symmetry planes is to calculate the 
value of the objective function for the trial estimate fa of the azimuth of a mirror 
symmetry plane. It would then be tested whether the value of the objective function for 
fa was zero, or acceptably close to zero; if it was, this would indicate that the azimuth of 
a mirror symmetry plane was, or was very close to, fa. If, however, the value of the 
objective function for fa was not zero or acceptably close to zero, then any suitable 
optimisation technique can be used to determine the azimuth of a mirror symmetry 
plane that minimises the objective function. 
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It should be noted that it is not necessary to investigate the objective function over the 
full 360° azimuth range. In the case of two orthogonal planes of mirror symmetry, any 
azimuth range of 90° should contain an azimuth value that minimises the objective 
function x 2 . In principle, therefore, investigating the objective function over an azimuth 
range of at least 90° should identify a mirror symmetry plane. Investigating the 
objective function over an azimuth range of less than 90° is not certain to identify a 
mirror symmetry plane, although it could do. However, if there exists a priori 
information about the likely orientation of a minor symmetry planes it would be 
possible to investigate the objective function over a small azimuth range around the trial 
estimated azimuth of a mirror symmetry plane. Such a priori information about the 
azimuth of a mirror symmetry plane might be derived, for example, from knowledge of 
the regional stress regime in a case where the azimuthal anisotropy results from the 
presence of fractures within the earth. 

The seismic data will have been acquired with a particular direction, for example north, 
defined as the azimuth zero. Once the azimuths of the mirror symmetry planes have 
been detected, it is possible to rotate the seismic data so that the azimuth zero is 
coincident with a mirror symmetry plane. This facilities subsequent processing of the 
seismic data, so allowing an improvement in the seismic image derived from the 
seismic data. 

Rotating the seismic data also allows extraction of further seismic attributes to 
characterise the seismic anisotropy. This may be done by a three-stage process in which 
initially a determination of the orientation of the symmetry planes is made using one 
azimuthally-dependent attribute of the seismic data. The data is then rotated so that the 
azimuth zero is coincident with a mirror symmetry plane. Finally, other attributes are 
extracted from the rotated seismic data. 

The invention has so far been described for data acquired using a surface seismic 
surveying arrangement of the type shown in Figure 1, in which the seismic source S and 
the receiver R are both disposed on the earth's surface. The invention is not limited to 
processing seismic data obtained by a surface seismic surveying arrangement, however, 
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and can be applied to all seismic data that show azimuthal anisotropy, regardless of how 
the data were obtained. For example, seismic data obtained using a VSP seismic 
surveying arrangement, in which the receiver R is disposed within a borehole, may be 
processed by a method of the invention, provided that the azimuth coverage of the 
acquired seismic data is sufficient. 

Figure 6(a) is a schematic plan view of a VSP surveying arrangement suitable for 
acquiring seismic data for processing by a method of the invention. A receiver R is 
disposed within a borehole B, and a plurality of source positions Si are defined around 
the receiver R on a circle centred on the borehole. In Figure 6(a),« a source position Si is 
defined every 10° of azimuth, but the azimuthal interval between adjacent source 
positions could be greater or less than 10°. 

In principle a separate seismic source could be located at each source location Si. It is 
generally more efficient, however, to use a "walkaround" technique in which a seismic 
source is disposed at one of the source locations Si and is actuated, and seismic data are 
acquired for that particular source-receiver azimuth. The source is then moved to 
another of the source locations Si and data are acquired at the new source-receiver 
azimuth. The process of moving the seismic source and acquiring data at the new 
source-receiver azimuth is repeated until data have been acquired for all desired values 
of the source-receiver azimuth. 

The invention can also be applied to reverse VSP seismic data, which is acquired using 
a seismic source disposed within a borehole and a receiver located on the earth's 
surface. Reverse VSP data suitable for processing by the method of the invention could 
be obtained using a seismic surveying arrangement corresponding to that of Figure 6(a), 
but with the source and receiver positions exchanged. That is to say, a seismic source S 
would be disposed within the borehole B, and a plurality of receiver positions Rj are 
defined around the source S on a circle centred on the borehole. 

Figure 6(b) illustrates another seismic surveying arrangement suitable for acquiring data 
for a number of source receiver azimuths. This is an Ocean Bottom Cable seismic 
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surveying arrangement, in which a seismic source array is moved through the water by a 
towing vessel. Seismic energy emitted by the seismic source array is detected by a 
static array of receivers disposed on the sea-bed. In the example shown in Figure 6(b) 
the array of seismic sources is a linear array, and the receivers are also arranged in a 
linear array. The source line is perpendicular to the receiver line, and this produces the 
common offset ring gather illustrated in Figure 6(b). The positions in Figure 6(b) where 
a thin line intersects the shot line or the receiver line indicate a source position or a 
receiver position respectively. 

In the example described above, the azimuth-dependence was determined for only one 
attribute of the seismic data. In general, however, seismic data acquired in a seismic 
survey may well contain more than one attribute that is azimuth-dependent. In this case 
two or more determinations of the orientation of the mirror symmetry planes may be 
made, each determination being made using a different attribute of the seismic data. 
The results of the separate determinations can be combined, for example using an 
averaging process, to obtain a final estimate of the azimuth of the mirror symmetry 
planes. 

Alternatively, two or more azimuth-dependent attributes may be combined, and the 
orientation of the symmetry planes may be determined from the combined attribute. 

The preferred embodiments described above relate to a case in which two orthogonal; 
symmetry planes exist within the earth. The invention is not limited to this case, 
however, and may be applied to the case where only a single symmetry plane exists. 
The existence of a single mirror symmetry plane requires that: 

= + (3) 

where / ( jf) is the value of the attribute at an azimuth <j>, and the symmetry plane lies at 
an azimuth In place of the objective function E given above in equation (2), a 
suitable objective function for the case of a single mirror symmetry plane is given by: 
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£'=t(/.-/W^o]) 2 



The objective functions E and E' given in Equations (2) and (4) above measure the 
square of the difference between the actual value of an attribute and the value predicted 
for that attribute using a trial estimate of the azimuth of a symmetry plane. Thus, the 
process of determining the azimuth of a symmetry plane that minimises the objective 
function E or E 1 is a "least squares" process. The invention is not limited to objective 
functions that measure the square of the difference between the actual value of an 
attribute and the value predicted for that attribute, however, and any objective function 
that provides a measure of the difference between the actual value of an attribute and 
the value predicted for that attribute can be used. 

As an example, an objective function that measures the absolute difference between the 
actual value of an attribute and the value predicted for that attribute could be used in 
place of the objective functions E or E\ Finding the azimuth of the symmetry plane 
that produces the "least absolute value" of the difference can be more reliable than a 
"least squares" process, particularly if the data contain noise. 

Thus, in the case of two orthogonal symmetry planes, a suitable "least absolute value" 
objective function would be: 



In principle, any objective function of the following general form could be used: 



(5) 



En = 2 I actual attribute - estimated attribute I N 



(6) 



where N is a positive number. 
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The invention has been described above with particular reference to seismic data. The 
invention is not, however, limited to seismic data but may be applied to any geophysical 
data that exhibits azimuthal asymmetry. For example, it could also be applied to 
borehole imaging logs constructed from, for example, acoustic and/or resistivity data. 

Figure 8 is a schematic block diagram of an apparatus 5 that is able to perform a method 
according to the present invention. 

The apparatus 5 comprises a programmable dataprocessor 6 with a program memory 7, 
for instance in the form of a read only memory (ROM), storing a program for 
controlling the data processor 6 to process seismic data by a method of the invention. 
The apparatus further comprises non-volatile read/write memory 8 for storing, for 
example, any data which must be retained in the absence of a power supply. A 
"working" or "scratch pad" memory for the data processor is provided by a random 
access memory RAM 9. An input device 10 is provided, for instance for receiving user 
commands and data. One or more output devices 1 1 are provided, for instance, for 
displaying information relating to the progress and result of the processing. The output 
device($) may be, for example, a printer, a visual display unit, or an output memory. 

Sets of seismic data for processing may be supplied via the input device 10 or may 
optionally be provided by a machine-readable data store 12. 

The results of the processing may be output via the output device 17 or may be stored. 

The program for operating the system and for performing the method described 
hereinbefore is stored in the program memory 7, which may be embodied as a 
semiconductor memory, for instance of the well known ROM type. However, the 
program may well be stored in any other suitable storage medium, such as a magnetic 
data carrier 7a (such as a "floppy disk") or a CD-ROM 7b. 
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CLAIMS: 

1 . A method of processing geophysical data, the data comprising at least a first set 
of geophysical data acquired with a first source-receiver azimuth and a second set of 
geophysical data acquired with a second source-receiver azimuth different from the first 
source-receiver azimuth* the method comprising; 

a) selecting an attribute of the geophysical data that is dependent on the source- 
receiver azimuth; and 

b) determining a value of the source-receiver azimuth about which the selected 
attribute exhibits mirror symmetry. 

2. A method of determining an azimufhal symmetry plane within the earth from 
geophysical data comprising at least a first set of geophysical data acquired with a first 
source-receiver azimuth and a second set of geophysical data acquired with a second 
source-receiver azimuth different from the first source-receiver azimuth, the method 
comprising the steps of: 

a) selecting an attribute of the geophysical data that is dependent on the source- 
receiver azimuth; and 

b) determining a value of the source-receiver azimuth about which the selected 
attribute exhibits miiror symmetry thereby to locate a mirror symmetry plane within the 
earth. 

3. A method as claimed in claim 1 or 2 wherein step (b) comprises determining a 
respective value of the selected attribute from each set of geophysical data. 

4. A method as claimed in claim 3 and further comprising the step of interpolating 
between the value of the selected attribute determined from the first set of geophysical 
data and the value of the selected attribute determined from the second set of 
geophysical data thereby to estimate the value of the selected attribute at a source- 
receiver azimuth intermediate the first and second source-receiver azimuths. 
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5. A method as claimed in claim 1 , 2, 3 or 4 wherein the determining step 
comprises 

c) calculating a function indicative of the variation between the value of the 
selected attribute at a first source-receiver azimuth (<f>- #?) and the value of the selected 
attribute at a second source-receiver azimuth - 0), where is the azimuth of a 
mirror symmetry plane, and 

d) determining a value of that minimises the value of the function. 

6. A method as claimed in claim 5 wherein step (d) comprises determining a value 
of fa that gives an value for the function that is lower than a pre-set value. 

7. A method as claimed in claim 5 or 6 wherein the function is: 

£ = S(//-/l^ + ^o]) 2 ^+a-/[180-^ + ^ 0 ]) 2 +(/;.-/[180 + ^-^]) 2 

M 

where/i is the value of the selected attribute at a source-receiver azimuth fa and /[a] is 
the value of the attribute at a source-receiver azimuth a. 

8. A method as claimed in any preceding claim and further comprising selecting 
another attribute of the geophysical data that is dependent on the source-receiver 
azimuth; and determining a value of the source-receiver azimuth about which the 
another selected attribute exhibits mirror symmetry. 

9. A method as claimed in claim 8 and further comprising the step of determining a 
value of the azimuth of a mirror symmetry plane within the earth from the value of the 
source-receiver azimuth about which the selected attribute exhibits minor symmetry 
and from the value of the source-receiver azimuth about which the another selected 
attribute exhibits minor symmetry. 

10. A method as claimed in any preceding claim wherein the or each attribute is 
selected from the arrival time, the polarisation angle, the transverse component, the 
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polarity, the temporal frequency, the amplitude or the linearity of acquired seismic 
energy. 

11. A method as claimed in any preceding claim and comprising the further step of 
controlling subsequent processing steps on the basis of the determined miiror symmetry 
plane. 

12. A method as claimed in any preceding claim and further comprising an initial 
step of acquiring at least a first set of geophysical data with a first source-receiver 
azimuth and a second set of geophysical data with a second source-receiver azimuth 
different from the first source-receiver azimuth. 

13. A method as claimed in any preceding claim wherein the geophysical data are 
seismic data 

14. An apparatus for processing geophysical data comprising at least a first set of 
geophysical data acquired with a first source-receiver azimuth and a second set of 
geophysical data acquired with a second source-receiver azimuth different from the first 
source-receiver azimuth, the apparatus comprising; 

a) means for selecting an attribute of the geophysical data that is dependent on the 
source-receiver azimuth; and 

b) means for determining a value of the source-receiver azimuth about which the 
selected attribute exhibits mirror symmetry. 

15. An apparatus for determining an azimuthal symmetry plane within the earth 
from geophysical data comprising at least a first set of geophysical data acquired with a 
first source-receiver azimuth and a second set of geophysical data acquired with a 
second source-receiver azimuth different from the first source-receiver azimuth, the 
apparatus comprising: 

a) means for selecting an attribute of the geophysical data that is dependent on the 
source-receiver azimuth; and 
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b) means for determining a value of the source-receiver azimuth about which the 
selected attribute exhibits mirror symmetry thereby to locate a mirror symmetry plane 
within the earth. 

16. An apparatus for performing a method as defined in any of claims 1 to 13. 

17. An apparatus as claimed in claim 14, 15 or 16 and comprising a programmable 
data processor. 

18. A storage medium comprising a program for a data processor of an apparatus as 
defined in claim 17. 
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